Shearing active gels close to the isotropic-nematic transition 
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We study numerically the rheological properties of a slab of active gel close to the isotropic-nematic 
transition. The flow behavior shows strong dependence on sample size, boundary conditions, and 
on the bulk constitutive curve, which, on entering the nematic phase, acquires an activity-induced 
discontinuity at the origin. The precursor of this within the metastable isotropic phase for contractile 
systems (e.g., actomyosin gels) gives a viscosity divergence; its counterpart for extensile (e.g., B. 
subtilis) suspensions admits instead a shear-banded flow with zero apparent viscosity. 

PACS numbers: 87.10.-e, 47.50.-d, 83.60.Fg, 47.63.Gd 



Active gels, including cell extracts and bacterial sus- 
pensions, are important and fascinating complex fluids 
HL 0, 0, 0, 0, 0, 0] • These systems contain self-propelled 
subunits which, on symmetry grounds, have generic or- 
dering tendencies that favor an isotropic-nematic (I-N) 
transition. Unlike conventional nematogcns, they are 
driven by a continuous energy flux and remain out of 
thermodynamic equilibrium even in the steady state. Ac- 
tive elements contribute to the stress as force dipoles that 
can be characterised as 'extensile' or 'contractile'. Rod- 
like extensile particles induce a dipolar flow away from 
the ends of the rods round to their long sides, whereas 
for contractile rods the flow direction is reversed. This 
means that an extensile, flow-aligning nematogen placed 
in a shear flow produces a flow that enhances the ap- 
plied shear; its contractile counterpart opposes the ap- 
plied shearing. 

In both cases, the quiescent bulk ordered phase (N) is 
unstable, but restabilizes at high shear rates Confin- 
ing walls can also restore stability 0, 0] , creating a finite 
threshold of activity for onset of spontaneous flow. A 
recent analysis 0] predicts moreover that the zero-shear 
viscosity of a contractile solution should diverge at the I- 
N transition. Thereafter, within the N phase, one expects 
a positive yield stress whose counterpart for an extensile 
nematic is, however, formally negative [l|. 

In this Letter, we study numerically the steady-state 
behavior of sheared active gels near the transition to a 
nematic phase. We discover a shear-banding instability 
in the extensile case which replaces a negative-viscosity 
zone by a window of 'superfluidity', accommodating a 
range of macroscopic shear rates at zero stress. We also 
find for contractile systems a viscosity divergence Q, at 
the (spinodal) threshold of the I-N transition. Near this, 
we find discontinuous shear thickening at a shear-induced 
I-N transition: though present in passive nematogens [1] , 
the latter is strongly perturbed by activity. The slab's di- 
mensions and boundary conditions can strongly influence 
this strikingly rich range of phenomena. 



We hope that, as well as offering insights into the 
highly nonlinear physics of active gels, our numerics will 
prompt new rheological experiments on soft active mat- 
ter such as bacterial suspensions and actomyosin gels. 

To gain our results we used a hybrid lattice Boltzmann 
(HLB) algorithm Q, whose approximations are shared by 
the hydrodynamic equations of motion (detailed below) 
of active fluids [3,0,0]; and also a finite difference (FD) 
algorithm that additionally neglects inertia [a, [9(. We 
assume translational invariance along flow and vorticity 
directions, reducing the 3D problem to ID (Io| . 

Governing Equations: In a coarse grained approach 
the local properties of an (apolar) active fluid can be de- 
scribed by a tensor order parameter, Q a p, whose largest 
eigenvalue, 2g/3, and its associated eigenvector give the 
magnitude and direction of the local oricntational order. 
The equilibrium physics of the passive system is then re- 
produced by a Landau-de Gennes free energy density T: 

(! ~ |- )~ljr ~ -^QapQpyQ-fa + ^(Qlp) 2 + y i^Qapfi 1 ) 

where tp controls the magnitude of the ordering. There 
is a first order transition between the N phase, sta- 
ble for (p > 2.7 and the I phase, stable for p < 2.7. 
At ip = 3.0 lies the spinodal point, the limit of the 
metastability of the I phase. In the last term of (Q]) 
K is an elastic constant (in the one constant approxi- 
mation) The equation of motion for Q a p is [l2j 
DtQafi = with Dt a material derivative describing 
advection by the fluid velocity u a , and rotation/stretch 
by flow gradients (see 0, E3|)- The molecular field is 
H a = — s ^ +((5 a /j/3)Tr ^ - and T is a collective rota- 
tional diffusivity. The fluid velocity obeys the continuity 
equation, and the Navier-Stokes equation 

p(d t + updp)u a = dp(TL a p) + (d a u/3 + d p u a ) (2) 

with p the fluid density, 77 a Newtonian viscosity, and 

n Q/ 3 = — Po5 a p + 2£(Q a p + -S a p)Q ie H je 
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The last term of ^ stems from activity, with C > 
(extensile) or ( < (contractile) [j], 0, Si ■ The pa- 
rameter £ depends on the geometry of active elements; it 
determines whether the material is flow-aligning or flow- 
tumbling. We restrict ourselves to the former by taking 
£ = 0.7. Our equations describe for simplicity apolar ac- 
tive gels, whereas many such gels (e.g., actomysin) are 
polar; however, we expect similar rheology, as the hydro- 
dynamic active stress has the same form in both cases [l| . 
Rheological experiments might also be possible on apolar 
active matter, such as melanocytes or vibrated rods 

Constitutive Curves: We first compute [IH curves (7(7) 
for the shear stress a at uniform imposed shear rate 
7 = d y u x . FigQ] shows this curve (and its stress con- 
tributions) for an extensile nematic. With our parame- 
ters, the passive terms in Ti xy are much smaller than the 
viscous term 777, so a ~ 777 — CQxyd)- I n both contrac- 
tile and extensile ncmatics Q xy has an upward disconti- 
nuity at the origin (see FigfT]); for contractile materials 
(C < 0) this creates an upward step in (7(7), or yield 
stress p|. The constitutive curve for the extensile case 
(C > 0, Fig[T]) instead has a downward step, or 'neg- 
ative yield stress', causing a zone of negative viscosity 
(cr/7 < 0) for I7I < 7*, which can be circumvented by 
splitting into shear bands of 7 = ±7*. The resulting com- 
posite flow curve 17(7), with 7 now the applied (mean) 
shear rate, connects 7 = ±7* with a horizontal 'tie line' 
at zero stress. The stability of such bands is far from cer- 
tain, especially if translational invariance is not assumed. 
However, our FD work (data not shown) confirms the 
existence of such states in the ID bulk (large L) limit. 
Figs. 5-7 of on confined slabs also offer evidence for 
this N/N banding scenario which, we argue below, should 
also control flow of extensile systems initialized in the I 
phase close the FN spinodal. 

Slab Geometry: We next consider a slab of active 
gel between parallel plates at y-separation L, with x- 
velocities 1L7L/2. Most of our HLB simulations were 
performed with L = fOO, T = 0.33775, K = 0.04, and 
C = ±0.005. These parameter values, and our results, are 
reported in LB units; such units map e.g. to Ax = 0.05 
//m and At = 0.67 /is with T = f Poise, n = 1.1 Poise, 
K = IQpN [ll| and C = 500 Pa. The only length scale in 
Eqs.CDE]) is then I = K 1 / 2 ~ 5nm. (This sets the scale 
for band interfaces, although these can be broader close 
to any second-order transitions.) Our HLB parameter 
represent molecular nematogens in a microfluidic-device- 
scalc slab a few /im across Q although very different 
values might arise in some biological contexts. 

Below we compare 'fixed' (anchored) boundary condi- 
tions, in which Q a p is specified on the walls, and 'free' 
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FIG. 1: For an extensile nematic (p — 3.0, £ = 0.7, n = 
0.567, £ = 0.1), as functions of 7 in steady state :(a) Passive 
order-parameter stress H^y = H xy + CQxy', (b) Total passive 
stress Yi^y = Hxy + 7 77; (c) Order parameter Q xy ; (d) Total 
shear stress (7(7) (dotted curve: the same for tp = 2.95.) 



ones which effectively set d y Q a p = 0, at y = 0, L. For 
shear-banding problems, the latter choice promotes con- 
vergence to bulk behavior our FD work uses this, 
and to the same end, chooses small £/L = 5 x 10~ 4 . The 
boundary condition on u a is taken as no-slip. 

Linear Contractile Rheology: We now report apparent 
viscosities 77ap P (7, L, ...) = {U xy +r]dyU x )/j. For homoge- 
neous flows, these give the bulk flow curve 77 app = cr( A f)/ A / 
and reduce to the zero-shear viscosity as 7 — > 0. 

We start with the zero-shear viscosity of a contractile 
gel in the metastable I phase, 2.7 < ip < 3.0. We observe 
a viscosity divergence at the latter (spinodal) point, as 
predicted in Q which, however, addressed strictly 2D 
nematogens for which binodal and spinodal coincide (as 
there is then no cubic term in ([T])). This is shown in 
Fig. [2^i, where the divergence is seen with free bound- 
ary conditions, for which the flow remains homogeneous, 
but is suppressed if Q a p is anchored at the walls with 
director along the flow (creating nonuniform shear near 
the walls) . In the latter case, 77 app increases linearly with 
L (data not shown; we have simulated up to L = 400) 
until it saturates, for ip < 3, at the free-boundary value. 
To understand this, note that surface anchoring lifts the 
rotational degeneracy of the N phase, so that deviations 
of the director n away from x feel a finite restoring force 
(that vanishes only as L — > 00). Hence the active shear 
stress, ~ —(qn x n y , remains of order 7 and there can be 
no yield stress, nor any viscosity divergence. 

Extensile Rheology: Consider next systems initialized 
in the metastable I phase of an extensile material (C = 
0.005). These have strikingly different behavior, shown 
in Fig. [2Jd. Not only does r7 app decrease strongly as the 
spinodal is approached: it becomes strictly zero in an in- 
terval </9 c (C, L) < <p < 3. We argue that this 'supcrfluid' 
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FIG. 2: (a) Apparent zero shear viscosity for an isotropic 
contractile gel in a slab under shear (solid and dashed lines). 
Solid curve(s): free boundary conditions; HLB data are indis- 
tinguishable from those derived from the bulk flow curve pl| . 
Dashed curve: fixed boundaries (director along flow), (b) Ap- 
parent viscosity in a linear regime for an active extensile gel 
in a slab under shear. Solid and dashed curves correspond 
to systems with free and fixed boundaries respectively. The 
insets show plots of Q xy and u x (solid and dashed lines) in 
the 77a PP = phase. 



window replaces one of negative viscosity on the constitu- 
tive curve, arising for ip > ip* (£) (say) , and similar to that 
already discussed for extensile nematic phases (Fig. [T|). 
To understand its continued presence here, note that as 
£p — > 3~ and at 7 = 0, dQ xy /d^ — ► —00; hence at small 

7 1 the flow curve has large negative slope (precursor to 
the step discontinuity in the nematic). As 7I increases, 
nonmonotonic features can arise, but the flow curve then 
joins that of the nematic, before stability is restored at 

7I = 7*. Close enough to the spinodal, therefore, the 
flow curve of the metastable I phase still has all the fea- 
tures needed for formation of bulk N/N shear bands. We 
then expect </? c (C, L — > 00) = ip*{Q] this is numerically 
consistent with our HLB and FD data. 

However, bulk shear bands do not provide a complete 
picture for the confined slab of Fig. [2b, where Q xy and u x 
remain relatively smooth on the scale of L (inset). With 
fixed boundary conditions, the zero-viscosity window is 
narrower and can disappear at small L. A complemen- 
tary interpretation of our N/N-bandcd state is as a 'spon- 
taneous flow' phase, predicted inside the nematic region 
of confined extensile materials 0, 0] ■ Our results suggest 
that, on increasing ip, a direct transition from metastable 
isotropic to the spontaneous flow phase is possible. 

Nonlinear Contractile Rheology: We again focus on the 
metastable I regime, 2.7 < (p < 3.0, whose passive coun- 
terpart shows shear thinning, and a shear-induced I-N 
transition [13 . In Fig. rj we plot HLB results for r7 app (7) 
in a confined slab of a contractile gel, with each shear 
rate 7 initiated from rest. We find a rich phenomenol- 
ogy, with strong dependences on boundary and initial 
conditions. Consider first the solid line, for ip = 2.99, 
with an initially isotropic state, and fixed Q al 3 = on 
the walls. For small 7 there is a linear regime of constant 
viscosity; the material remains in the I phase. There is 
then a sharp increase in ?7 app , i-e., strong shear thickening 
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FIG. 3: Apparent viscosity vs. 7 for ip = 2.99 (solid, dashed 
and dotted line) and ip = 2.8 (dot-dashed line) for contractile 
gels under shear. Anchored BCs are considered, with Q a p — 
(solid, dot-dashed) or with n || x (long-dashed and dotted 
lines). For the dotted line, the bulk was also initialised with 
nematic order. 



(discontinuous at large L/£), with onset of nematic order 
oriented near but above the Leslie angle 8 — 9*, and a 
drastically non- Newtonian flow profile (Fig. [3l inset). At 
larger 7, 9 decreases and 77 app falls towards 77. The dashed 
curve shows the same system but with fixed nematic or- 
dering and n || x at the walls. For small 7 the shear 
flow is now mainly confined to wall layers; the director 
has 9 ~ throughout the slab. In the wall layers, the 
backflow term Ti a p acts to restore the ordering, creating 
extra dissipation and enhanced 77 app at small 7. (This 
scenario is reminiscent of a passive nematogen flowing 
in the permeation mode close to the isotropic-cholesteric 
transition [3], where energy from the applied shear is 
converted into a flow that opposes the shear and reduces 
the order.) Shear thickening is now accompanied by an 
abrupt director transition, near the slab centre. As 7 
is raised, 9 decreases, ?7 app falls, and a Newtonian flow 
profile is restored. If, instead, the system is initialised 
in the nematic phase (dotted curve) the jump in 9 dis- 
appears, as does the peak in 77 app . Finally, we consider 
an initial I state further from the spinodal at ip = 2.8, 
with Q a p = at the wall (dot-dashed line). This system 
remains isotropic for longer: ?7 app is smaller in the linear 
regime, remains there to higher 7, and shows a smaller 
peak at the shear-induced I-N transition. 

These results for confined slabs are again illuminated 
by computing the homogeneous bulk flow curve o-(j) [l3j |. 
For C < this has a low viscosity I branch, terminating 
at some finite flow rate, and a high viscosity, shear thin- 
ning N branch (reflecting the presence of a yield stress). 
Fig.|4]shows both the flow curves (inset) and the viscosity 
77(7). Discontinuous shear-thickening arises on tracking 
up to the end of the I branch, then jumping vertically 
to N. (On decreasing shear rate, N can be tracked down- 
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FIG. 4: Viscosity and (inset) flow curves for active con- 
tractile gels in bulk homogeneous shear: <p = 2.8, 2.9, 2.99 
(dot-dashed, dashed, solid). Vertical line segments signify 
the end of the I branch and a jump up to N. Data with 
nematic anchoring as in Fig. [3] is shown for an initially 
isotropic (squares) or nematic (circles) bulk. 
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FIG. 5: Dynamic phase diagram in the (C,wo) plane (with 
Vo — 7-L the wall velocity; L = 100) for an active liquid crys- 
tal close to the isotropic-nematic transition (<p = 2.695, r\ = 
0.28). The phases are labelled as follows: N= nematic, I = 
isotropic, Bl = I/N banded, B2 = N/N banded. 



kincsin- microtubule solutions. 



wards to low 7 and likely remains metastable.) Once 
again, there are features of the confined behavior, such 
as boundary-condition sensitivity, that this bulk picture 
cannot explain. 

We finally consider the effect of activity on the I/N 
shear banding present already in passive nematogens 
8]. To this end, we fix tp = 2.695 and K = 0.01; in the 
passive system (( = 0) there is a range of shear rates 7 
for which the fluid splits into I and N bands of high and 
low viscosity. Activity has a significant impact: as shown 
in Fig. [51 it widens or narrows the banding window 
according to the sign of £. (For larger <j> the bands pinch 
off giving a direct I/N transition, as reported above, in 
sufficiently contractile materials.) The trend concurs 
with expectations from the bulk flow curve: (7(7) for the 
N phase is lowered by extensile and raised by contractile 
activity, allowing the former system to band more easily. 

To conclude, we have solved the hydrodynamic 
equations of motion of an active nematogen close to 
the I-N transition subjected to a shear flow, assuming 
translational invariance along flow and vorticity direc- 
tions. For contractile gels and free boundary conditions 
there is a divergence of the apparent viscosity at the 
I-N spinodal Extensile gels, in contrast, enter a 

zero-viscosity phase of N/N shear bands as the spinodal 
is approached. In the nonlinear regime, both extensile 
and contractile materials show non-monotonic effective 
flow curves, with details strongly dependent on initial 
and boundary conditions. Finally, extensile (contractile) 
activity stabilizes (suppresses) I/N shear banding. We 
hope our results will spur experimental investigations of 
the rheology of active gels, such as extensile bacterial 
suspensions, or scmidilutc contractile actomyosin or 
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